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Abstract. We present the analysis of Rossi X-ray Timing Explorer (RXTE) observations of the turn-on phase of a 35 day cycle 
of the X-ray binary Her X- 1 . During the early phases of the turn-on, the energy spectrum is composed of X-rays scattered into 
the line of sight plus heavily absorbed X-rays. The energy spectra in the 3-17 keV range can be described by a partial covering 
model, where one of the components is influenced by photoelectric absorption and Thomson scattering in cold material plus 
an iron emission line at 6.5 keV. In this paper we show the evolution of spectral parameters as well as the evolution of the 
pulse profile during the turn-on. We describe this evolution using Monte Carlo simulations which self-consistently describe the 
evolution of the X-ray pulse profile and of the energy spectrum. 
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1. Introduction 

Her X-1 is one of the best understood X-ray binary sys- 
tems showing a variety of long and short term periodicities. 
The X-ray pulsar spins with a 1.24 s period and moves in a 
1.7 day almos t circular orbit around its companion HZ Her 
jXananbaum et al. 1972). Both effects cause a modulation of 
the observed flux in optical as well as in the X-rays. In addi- 
tion the X-ray light curve of Her X-1 shows a long term 35 
day intensity variation. This modulation is the best evidence 
for an inclined, processing, and warped accretion disk in an 
X-ray binary system. The origin of the warping is not yet fully 
understood, but may be caused either by r adiation driven accre- 
tion disk winds JSchan dl & Me veJI 19941) or by radiation pres- 
sure jMalonev & B egel'mar3 Il997i) . The processing motion of 
the disk can be understood in the context of tidal interaction 
and/or as a consequence of non van ishing torques acting on 
the disk, e.g. due to a coronal wind ( Schandl & Meve]"l994j 
ISchmdi et al,J997: Shakura etal. 1998; Kets aris et al.,.200li) . 
Because of the high inclination of the system, the disk period- 
ically blocks the line of sight to the neutron star during about 
60% of the 35 day cycle. 

The observed 35 day light curve shows t wo maxima in 
intensity: the "main-on" and the "short-on" jGiacconi et alJ 
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119731) . Following the often adopted baseline model of 
Her X-1 jKatz 1973: Schand l & Mever 1994 : Scottet^lv 2600I 
iKetsaris et alJl20o'l[lLeahvll2004l and references therein), the 
~ 12 day long main-on starts when the outer rim of the accre- 
tion disk opens the line of sight to the central neutron star (see 
Fig.|5}- Subsequently, at the end of the main-on the inner edge 
of the accretion disk covers the line of sight. The second maxi- 
mum in intensity occurs as soon as the inner edge of the accre- 
tion disk uncovers the line of sight during the progression of the 
35 day cycle. This phase is called short-on where only ~35% 
of the main-on intensity is measured. The states in between the 
short-on and the main-on are called "low-states " where the in- 
tensity drops to ~3% of the main-on intensity dScott & Leahvl 
fl999). The transitions from the low-states to the main-on and 
short-on are called "turn-on", while the decline in intensity at 
the end of the main-on and the short-on generally is called 
"turn-off". This periodic behavior is irregularly i nterrupted by 
anomalous l ow states as it was observed in 1983 (IParmar et al.l 
|1985'), 1993 /Vrtile k et alJl994tlVrtilek & Ch eng'l9 96|), 1999 
garmar et al. 1999: ICoburn et al.l l2000t [Oosterb roek et alJ 
l200ll) . and recently in 2004 iBovdetal] E004I) . During an 
anomalous low state the maximum X-ray flux typically drops 
to 1-3% of the main-on flux. 

While the end of the main-o n has been the subj ect of pre- 
vious Ginga observations (e.g., IPeeter et al.lll998i) . observa- 
tional data on the spectral evolution during the start of the 
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Fig. 1. From top to bottom; RXTE PCA and RXTE ASM count rates for the time of the turn-on. The beginning and the end of 
an echpse are marked by dashed Hnes. The numbers identify individual RXTE orbits starting with orbit 00 to orbit 31. 



main-on has been rar e. Observations of t u rn-ons of Her X- 1 
have b een p resented bvlBecker et alJ lll977l) . lDavison & FabianI 
( Il977l) . and lParmar et alJ(ll980l) with different instruments. All 
of these observations show a strong indication that the flux dur- 
ing the early stages of the turn-on is composed of heavily ab- 
sorbed plus scatted radiation. In 1997 September we observed 
a complete turn-on of the 35 day cycle of Her X-1 in a two day 
continuous observation with the Rossi X-ray Timing Explorer 
(RXTE). In this paper we present results from the spectral and 
temporal analysis of this observation. Further we describe these 
data with a physical model that can reproduce the spectral fea- 
tures and temporal evolution of the pulse shape seen in this 
observation. For our physical model we assume that the varia- 
tions of the observed spectrum are due to the varying column 
density caused by cool gas of the outer rim of the accretion disk 
and due to radiation scattered from (ionized) gas sandwiching 
the accretion disk (an accretion disk corona). 

The remainder of this paper is structured as follows; In 
Section|2we give a short description of our observation and the 
extraction of the data, before we describe our spectral model 
and the results of the spectral analysis in Section|3l we present 
the results of our analysis of the evolution of the pulse profile 
in Section^] In Sectionlsjwe introduce a method to determine 



the amount of absorbed and scattered radiation by simulating 
the influence of a scattering medium on the shape of the pulse 
profile using Monte Carlo simulations. We summarize this pa- 
per in Section |6l and propose a model of the outer accretion 
disk rim which can explain the spectral behavior as well as the 
evolution of the pulse profile. 

2. Observation and data reduction 

We observed Her X-1 for two days on 1997 September 13/14 
with RXTE. Fig. ^ shows the RXTE Proportional Counter 
Array (PCA) light curve of the entire observation, together 
with the light curve measured simultaneously by the RXTE 
All Sky Monitor (ASM). Note that during the entire observa- 
tion all five PCUs of the RXTE PCA were active and therefore 
throughout this paper the count rates are consistently given in 
cts/s for five PCUs. During the turn-on an eclipse took place 
around MJD 50705.5. Furthermore, two dips were detected; 
a pre-eclipse dip around MJD 50705.3 and an anomalos dip 
around MJD 50704.8. The gaps in the light curve are due to 
Earth occupations and SAA passages during individual RXTE 
orbits. The exposure times, mean count rates, and dates of ob- 
servation are given in Table[nfor each RXTE orbit. 
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Table 1. Observing log of the Turn-On observations of Her X- 
1. 



Obs. 


Date 
[MJD] 


[sec] 


[counts s"'] 


00 


50703.979 


3300 


62.1 ±0.2 


01 


50703.979 


3400 


73.3 ±0.3 


02 


50703.979 


3300 


100.0 ± 0.3 


03 


50703.979 


3300 


130.0 ± 0.3 


04 


50703.979 


3200 


195.4 ±0.3 


05 


50704.312 


3000 


280.8 ± 0.4 


06 


50704.312 


2600 


402.7 ± 0.5 


07 


50704.452 


2300 


469.2 ± 0.5 


08 


50704.452 


2100 


395.0 ± 0.5 


09 


50704.452 


1700 


501.7 ±0.6 


10 


50704.727 


1200 


158.8 ±0.5 


11 


50704.727 


300 


148.7 ± 1.0 


12 


50704.727 


2100 


160.6 ± 0.4 


13 


50704.727 


700 


156.7 ± 0.7 


14 


50704.860 


2100 


358.4 ± 0.5 


15 


50704.860 


3400 


813.6 ±0.5 


16 


50704.860 


3300 


1197.0 ±0.6 


17 


50705.045 


3400 


1298.1 ±0.7 


18 


50705.045 


3200 


1293.7 ± 0.7 


19 


50705.045 


3100 


325.6 ± 0.4 


20 


50705.312 


3000 


329.3 ± 0.4 


21 


50705.381 


2700 


34.9 ± 0.3 


22 


50705.591 


1800 


73.2 ± 0.4 


23 


50705.659 


1500 


1494.1 ± 1.1 


24 


50705.726 


1200 


1568.2 ± 1.2 






9 inn 


1 H\i ^ + n Q 


26 


50705.793 


3300 


1612.5 ±0.7 


27 


50705.793 


3400 


1697.5 ± 0.7 


28 


50705.793 


3400 


1699.1 ±0.7 


29 


50705.793 


3100 


1787.7 ± 0.8 


30 


50706.096 


3100 


1810.4 ±0.8 


31 


50707.313 


1900 


1828.2 ± 1.0 



3 

Table 2. Systematic eiTors applied to the PCA data to account 
for uncertainties in the PCA response matrix. 



Channel 


Energy Range 


Systematics 


0-15 


2-8 keV 


1.0% 


16-39 


8-18 keV 


0.5% 


40-57 


18-29keV 


2.0% 


58-128 


29-120 keV 


5.0% 



imum flux in the main pulse of the profile in the energy range 
13-19 keV. From each profile we subtracted the unpulsed flux 
and normaUzed the count rate to the maximum. Fig. |6l shows 
the evolution of the pulse profiles in different energy ranges 
over the time of the whole turn-on. We have omitted those or- 
bits during which the pulse profile shows no remarkable vari- 
ation compared to the previous or following orbit, i.e., the or- 
bits 01, 06-08, 16-18, and 25-30. 

For the spectral analysis, pulse phase averaged PCA and 
HEXTE spectra were extracted for each individual RXTE or- 
bit. To minimize the background, we have chosen good-time 
intervals (GT I) with an "electron -ratio" of all PCUs less than 
0.1 (see e.g.. lWilmse7aDll999l) . All spectra are background 
and dead-time corrected. The data of orbits 10-14 and 19-22 
were omitted for the analysis because our spectral model is not 
applicable during the times of the dips and the eclipse. We si- 
multaneously fitted our spectral model described in section lTTI 
to both the HEXTE and PCA data of each orbit. The system- 
atic uncertainties of the response matrix of the PCA assumed 
for the spectral analysis are given in Table|2 



3. Evolution of spectral parameters 

Before we present the results of our spectral fitting in Sect. 13. 2! 
we give an introduction to the complex spectral model used to 
describe the data. 



Exposure times shown are rounded to the closest 100 sec. The count Spectrsi model 

rate is background subtracted. j 1 

As earlier observations already have shown ([Davison & Fabian! 

Il977t IParmar et al.lll980t iBecker et al.lll977l) . a combination 

For the detailed analysis, we extracted PCA light curves of direct, scattered, and absorbed photons is observed during 

for all RXTE orbits listed in Table [T] with a time resolu- the turn-on. Therefore, we used a partial covering model which 

tion of 16 ms. Light curves were extracted for five energy combines both, scattering and absorption, to fit the data over the 

bands: 2.0-4.5 keV 4.5-6.5 keV 6.5-9 keV 9-13 keV and 13- of '^he turn-on. As components for the spectral model we 

19keV. After correcting the photon arrival times with respect used an exponentially cutoff^ power-law /poweiC^) • /highecutC^) 

to the solar systems barycenter and for the orbital motion of as implemented in XSPEC, a cyclotron line feature /cyc(£) at 

the neutron star, we determined the pulse period of Her X-1 39keV using the cyclabs model of XSPEC, and a Gaussian 

by folding the data of the 13-19keV energy band using a emission line Ip^{E) fixed at 6.4keV. For an analytic descrip- 

maximization test. The resulting pulse period is Pspin = tion of the individual model components we refer the reader to 

1 .2377291(2) s (MJ D 50708. 199), which is consi stent with ob- ^e XSPEC manual (|Arnaud & Dorman|2002|). 

servations of, e.g., iDalFiume et al.l ( Il998l) and ICoburn et al.l We combined these four spectral components to provide a 

1I2OOOI) . Subsequently we folded all light curves with this pulse basic continuum model, which is then observed through an ab- 

period to obtain a pulse profile for each energy band and RXTE sorber partially covering the continuum source. In an analytic 

orbit. Pulse phase O1.24 - was defined as the time of the max- form the final spectral model describing the photon spectrum 
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Fig. 2. Schematic illustration of the geometry during a turn-on 
of Her X-1 (not to scale). The accretion disk rim, the accre- 
tion disk corona, and the neutron star are shown. Three differ- 
ent components contributing to the overall observed flux are 
indicated: radiation absorbed by the cold material of the accre- 
tion disk (photons following beam A), radiation absorbed by 
the accretion disk rim and further reduced in flux by Thomson 
scattering (photons following beam B), and radiation scattered 
by the corona into the line of sight of the observer (photons 
following beam C). 

can then be written as 

HE) ^ B ■ (I + a(E)) ■ (/()(£) /Fe(£)) (1) 
Primary Spectrum 

with 

a(E) = . c ■ e_^^ (2) 

Thomson-Scattering Absorption 

and 

hiE) = /powei(£') ■ 4ighecut(£') ' ^cyc(£') (3) 

where crj is the Thomson cross-section and where the bound- 
free absorption c ross sections, cThf , a re those used in the tbabs 
model of XSPEC jWilms et alJ200(]l) . For the fitting the absorp- 
tion and scattering part in Eq. (|2} uses the same A^h- For the 
Thomson scattering part in Eq. (0 the electron density cal- 
culates from A^H - 1.21 A^ c assuming material of solar abun- 
dances dWilms et a l."2000^. The constant C in Eq. (|3 defines 
the relative ratio of absorbed and scattered radiation to the un- 
affected radiation. A larger value of C implies a larger degree 
of absorbed and scattered flux. During the remainder of this pa- 
per the unaffected model component will be called MC 1 and 
the component influenced by absorption and Thomson scat- 
tering will be referenced as MC 11. For the model compo- 
nent MC 1 we neglected Galactic absorption since the Galactic 
A^H = 1-79 10^" cm"^ is small compared to the A^h of model 
component MC 11. A schematic picture illustrating the geomet- 
ric situation during the turn-on is shown in Fig.|2] The positions 
of the accretion disk and an accretion disk corona relative to the 
neutron star are shown. Photons following beam A, marked by 
a dashed fine, are scattered in the corona and can be partially 



absorbed in the accretion disk rim (beam C, dash dotted line). 
While a certain fraction of photons are blocked by the accretion 
disk (beam B), the outer parts of the accretion disk are optically 
thin for X-rays. Therefore, the spectral distribution of photons 
following beam D will show strong signature of photoelectric 
absorption. To simplify matters, photons reaching the observer 
directly without being modified by either the accretion disk nor 
the accretion disk corona are not shown in Fig|2] Comparing 
this geometric interpretation with the spectral model given in 
Eq. (|2j allows the following interpretation: 

- The unaffected model component I{E) - Io{E) + iFe(E), 
called MC 1 from now on, accounts for photons following 
beam A and photons reaching the observer directly (this 
case is not shown in Fig.|2j. 

- The model component modified by photoelectric absorp- 
tion and Thomson scattering I{E) = a{E) ■ {Io{E) + Ipf.{E)), 
called MC 11 from now on, represents the spectral distribu- 
tion of photons following beam C and beam D. 

We emphasize that it is not possible to use a physically more 
realistic spectral model which treats scattered, absorbed, and 
direct flux separately. The problem lies within the nature of 
Thomson scattering: for photon energies E < lOkeV {E «: 
mgC^) and when the influence of multiple scatterings can be ne- 
glected (r < 5), scattering of photons by stationary and free 
electrons can be treated as elastic and consequently energy in- 
dependent (classical Thomson approximation). This makes it 
impossible to separate, e.g., the direct flux and the flux contri- 
bution from photons scattered into the line of sight (beam A) 
using spectral analysis alone. 

Using this simplified partial covering spectral model all 32 
phase averaged PCA and HEXTE spectra were fitted in the en- 
ergy range 2.9-18 keV (PCA) and 15-100 keV (HEXTE). In a 
first iteration we fitted the data with all parameters free except 
the power-law index a which was kept fixed at 1.068. This is 
an average value for a derived from the data with high count- 
ing statistics towards the end of the turn-on (orbits 15-31). 
The results for the remaining free fit parameters are listed in 
Table lA.f1 This analysis reveals that the folding energy iScutoff. 
Efoid, Ecyc, (T eye, Efg, and crpj. show no significant variation 
over the duration of the turn-on. Therefore, these values were 
kept fixed at their mean values (see Tabled for the further 
analysis. Leaving these values fixed allows us to determine the 
variation of the remaining free parameters for the time of turn- 
on, which are the neutral column density A^h, the ratio C, the 
normalization of the power law Apl, and of the iron emission 
line ALine- The results are shown in Fig.|3and the correspond- 
ing fit parameters are given in Table lA.21 

3.2. Spectral Parameters 

Assuming a simple geometrical model of a turn-on where the 
outer edge of the accretion disk opens the line of sight to the 
central neutron star as indicated in Fig. |2 the neutral column 
density A^h is expected to gradually decrease when the flux in- 
creases. As can be clearly seen from Fig.|3^, the behavior of ob- 
served A^H is contr ary to this sirnple rn odel: As we have already 
mentioned earlier dKuster et al.l 19991) . A^h increases during the 
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Fig. 3. From top to bottom: (a) A^h [lO^^cm^^], (b) ratio between the model component MC II and MC I as defined in Eq. Q, 
and (c) normalization of the iron line Aiine- The uncertainties are +lcr. 



Table 3. Parameters fixed to their mean value for the fitting of 
the spectra over the time of the turn-on. 



Parameter 


Value 


a 


1.068 


^cutoff 


21.5 keV 


-£fold 


14.1 keV 


F 


39.4 keV 


f7"cyc 


5.1 keV 


Epe 


6.45 keV 


0"Fc 


0.45 keV 



first four RXTE orbits, reaches a maximum during orbit 5, and 
then declines until orbit 17 where it becomes untraceable. This 
evolution goes in parallel to the degree of scattered and ab- 
sorbed radiation, C, shown in Fig. which increases until 
orbit 09 and afterwards starts to decline. The turning point in 
the progression of C matches almost exactly the maximum of 
A^H- In contrast to A^h and C, the normalization of the Fe emis- 
sion line, Aune reproduces the progression of the count rate and 
increases during the whole of the turn on (Fig. |3j;)- We will 
come back to the interpretation of this behavior of A^h and C in 
Sect.m 



4. Evolution of the pulse profile 



4.1. Pulse variation depending on disk phase 



Further insight into the physics of the turn-on comes from the 
substantial changes in shape and amplitude of the X-ray pulse 
with 35 day phas e. One of th e earliest studies of these changes 
was presented bv lBail ( Il98ll) . who observed that the hard cen- 
tral peak and the soft traihng shoulder of the Her X-1 pulse 
profile are affected differently over the time of the turn-off. He 
interpreted this effect by a time dependent covering of the two 
polar emission regions on the neutron stars surface by the in- 
ner accretion disk rim and its corona. Since the thickness of the 
inner accretion disk rim and the neutron star are of the same or- 
der of magnitude, they have similar angular size. Consequently, 
the observed flux from the two neutron star poles can be af- 
fected differently in time by the material of the inner disk rim. 
Building upon this work and on observations of th e change of 
the pu lse-profile during the end of the 35 day cvcle. lScott et al.l 
(I2OOOI) were able to present a refined geometric model explain- 
ing these changes in more detail. 
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1.0 

Pulse Phase 

Fig. 4. Nomenclature for the different features of the Her X- 1 
pulse profile. This pulse profile is taken from orbit 30 in the 
energy range of 9.0-13.0 keV. 



4.2. Observed pulse variation 

For the following discussion, we use the nomenclature given 
in Fig. |?]for the various features of the Her X-1 pulse pro- 
file. Fig. |5] shows how the relative strength of these features 
change with energy. The most distinct variation apparent, is 
the decreasing flux of the soft leading shoulder with increas- 
ing energy. This change leads to a double peaked structure 
consisting of the hard central peak and the soft leading shoul- 
der that clearly can be identified in the energy range 6-12 keV 
jPeeter et al.ll998h . 

During the early phases of the turn-on, strong photoelec- 
tric absorption and Thomson scattering will modify the pulse 
shape (Fig.|6}. Thomson scattering in a hot plasma causes the 
broadening of the pulse profile in all energy bands which leads 
to an almost sinusoidal pulse shape during early times of the 
turn-on (orbits 00-05). During later phases of the turn-on (or- 
bits 05-09) only the lowest energy channels are affected by 
strong noise. This implies that during this phase energy depen- 
dant photoelectric absorption is the dominant process. The be- 
havior during orbits 05-09 is very similar to the situation in or- 
bits 10-14 during which the anomalous dip took place, which is 
presumably caused by cold material located at the outer rim of 
the accretion disk, crossing the line of sight to the neutron star 
(ghakura et al. 1993P- During the eclipse phase (orbit 21-22) 
all energy channels are affected by strong noise and no pulsa- 
tion was detected. 

These effects of scattering and photoelectric absorption can 
also be seen in the behavior of the pulsed fraction with time 
(Fig-E}- It is clearly visible that F pulsed increases more rapidly 
in the high energy bands. At lower energies, the pulsed flux is 
suppressed towards the beginning of the turn-on, similar to the 
situation observed during egress of the anomalous dip, where 
the pulsed flux increases faster at higher energies. At the end 
of the turn-on, after orbit 23, the pulsed fraction is almost con- 
stant. Contrary to later phases of the main-on and the turn-off 
the intrinsic pulse shape does not change significantly over the 
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Fig. 5. Energy dependence of the pulse profile. The pulse pro- 
file of orbit 31 is shown in eight different energy bands. All 
profiles are normalized to unity. Each profile, except the profile 
of the energy range 2.0-3.6 keV, is shifted by 0.5 in the y direc- 
tion relative to the previous profile. Note, the energy dependant 
change of the relative intensity of the soft leading shoulder to 
the hard central peak, which results in a double peaked struc- 
ture close to pulse phase 1 .0. This feature is most pronounced 
at energies between 10-14keV 



time of the tu rn-on. This result is in agre e ment with earlier find- 
of, e.g., 'Gruber et al.' ('1980"), 'Bal ( Il98ll) . iTriimper et all 



ings 



e.g., 

(11986ft . or iDeeter et al. ( 199&) . The pulse profile observed at 
the beginning of the turn-on can be interpreted, therefore, as a 
main-on pulse profile modified by the influence of photoelec- 
tric absorption and scattering. 



5. Simulating pulse variation 

To quantify the effects of scattering and photoelectric absorp- 
tion on the change of the X-ray spectrum and the pulse profile 
shown in Sects. |3and|2 we now turn to Monte Carlo simu- 
lations of the radiation transport of the pulsar's flux in a scat- 
tering medium. For the simulations we assume that the pulse 
shape seen at the end of the turn-on (orbit 30), is the intrin- 
sic pulse shape caused by the emission characteristic of the 
neutron star and is not changing over the time of the turn-on. 
Furthermore, we assume that the smearing of the pulse profile 
and the change in spectral shape at the beginning of the turn- 
on are solely caused by scattering and photoelectric absorption 
in the medium covering the line of sight to the neutron star 
With these assumptions we can use the pulse profile observed 
in orbit 30 as a "template" profile and investigate the effects of 
a scattering and absorbing corona on the pulse shape depend- 
ing on A^H and the size of the scattering region. This is done 
via Monte Carlo simulations, as described in the following sec- 
tions. 
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Fig. 6. Evolution of the pulse profile as observed over the time of the turn-on. All profiles are normalized to unity at the maximum 
of the main pulse, after subtraction of the off-pulse constant flux. Pulse phase is defined as the maximum of the main pulse in 
the energy band of 13-19keV. 
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Fig. 7. Pulsed fraction over the time of the turn-on in different energy ranges. To determine the pulsed fraction the pulse profiles 
have to be smoothed to reduce Poisson noise. The pulsed fraction is given as percentage of the pulsed flux compared to the non 
pulsed flux. 



5.1. Monte Carlo Simulations 

Assuming a point source emitting the intensity /(?o) at time to, 
the intensity at infinity observed at time t can be written as 

r{t)= f G{t,to)I{to)dto (4) 

where G(f, fo) is the scattering Green's function, i.e., the appro- 
priately normalized solution of the time-dependent equation of 
radiation transfer through the scattering and absorbing medium 
for a 5-function pulse of light emitted at time fo. 

First analytical approaches to determine the Green's 
function for scattering in a cold coron a were based 
on the fundamental method dev e loped bv iLightman et alJ 
lll98ll). iBrainerd & LambI lll987h . iKvlafis & KlimisI l ll987h . 
and iKvlafis & PhinnevI ( 1 19891) found analytical solutions for 
G(f, fo) for simple geometries, such as the diffusion of photons 
out of a spherical shell surrounding a central point source. For 
the general case, analytical results are difficult to obtain and 
one has to resort to numerical solutions instead. Here, we use 
a modification of t he linear Mo nte Carlo code based on the 
method of weights (ISoboil99l|) that we have previously used 
to compute G(f, fo) for the case of a hot Comptonizing plasma 



( lNowaketalJll999l) . In our simulations we consider Compton 
scattering (using the Klein-Nishina cross section), photoelec- 
tric absorpt ion from material of solar abundance using the cross 
sections of IVerner et alJ ( Il996l) . and fluorescent line emi ssion 
(using the fluorescence yields of iKaastra & Mew'3ll993h . We 
model the propagation of photons through a plane-parallel slab 
with thickness d and hydrogen column A^h (corresponding to a 
certain optical depth for electron scattering, Tes) that is illumi- 
nated by a source at infinity. We consider both, neutral and fully 
ionized slabs. Output of the simulation is G(f, fo) as a func- 
tion of A^H, d, and energy band and the angle-dependent photon 
spectrum of photons leaving the slab. We normalize the time to 
the light crossing time of the slab. 

As an example, Fig.|8lshows G(f, fo) for a neutral slab with 
A^H = 1.010-2 cm-2, 1.25 1023 cm^^, and 2.75 lO^^ cm^^, 
while Fig.|5]shows the same for a fully ionized slab. For Tgs < 5 
the Green's function G(f, fo) in Fig.|S]and Fig.|5]is dominated 
by photons crossing the slab on a straight line without being 
scattered from electrons. These photons contribute to the peak 
apparent at diffusion times equal to the light crossing time of 
the slab (fjiff = 1) and only a small number of photons dif- 
fuse out of the slab after this initial peak. For increasing elec- 
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2 4 6 

Light crossing time [slab thickness] 

Fig. 8. G{t,to) in the energy range 1.0-70.0keV for a neu- 
tral corona and low optical depths of Tes - 0.01 (A^h = 



1.010^2 ^,^-2^^ ^ 0.14 (A^H = 1.25 10^" cm-^), and 
= 0.44 (A^H = 2.75 10^"^ cm^^) 'phe peak at f = 1 is caused 
by photons crossing the slab without scattering, the break at 
f = 2 is caused by photons scattering at most one time be- 
fore leaving the slab. Inset: Spectrum emerging from the slab 
for the same columns. Note the emergent fluorescent emission 
lines for the larger optical depths. 



tron optical depth, the number of diffusing photons increases 
significantly since the mean number of scatterings per photon 
increases approximately as ~tI^. As a consequence, the max- 
imum of G{t, to) moves towards higher diffusion times, un- 
til G{t, to) is dominated by photons scattered multiple times 
(Tes > 5). The energy loss per electron scattering event can 
be estimated for a 15 keV photon to AE x 0.44 keV. This im- 
plies that, as soon as the mean number of scattering per photon 
increases above ~ 15-20 (res > 4), photons of the energy range 
between 10-20keV are redistributed to the energy band of 2- 
lOkeV. This redistribution effect is the origin of the changing 
spectral shape with increasing optical depth visible in Fig. O 
(inset). As a consequence, the cut-off energy moves towards 
lower energies and a bump at energies between 1-5 keV arises. 
For even larger optical depths, Tes > 10, this bump slowly van- 
ishes and the flux above 6 keV decreases rapidly. The overall 
spectral shape is then given by the left most spectrum shown in 
the top panel of Fig.|51 

Considering a neutral medium, absorption plays the domi- 
nant role over the temporal effects of Compton scattering. For 
a neutral medium with Tes ^ 1, almost all flux below lOkeV is 
suppressed'. At such optical depths the effects caused by dif- 
fusion time are still negligible since the mean number of scat- 
terings per photon is close to unity. Thus, to achieve noticeable 
changes in beam shape a high fraction of ionized material is 
needed. On the other hand a high fraction of neutral material 



' For material with solar abundance, the photoionization cross sec- 
tion equals the Thomson cross section at ~10keV. Due to its E^^ pro- 
portionality the optical depth below this energy will therefore be al- 
ways higher than the electron optical depth. 



2 4 

Light crossing time [slab thickness] 

Fig. 9. G{t,to) for the same spectral distribution as shown in 
Fig Elbut for a fully ionized medium and higher optical depths. 



only weakly alters the beam shape but reduces the flux at low 
energies. 

Fi g .1 1 Oldemonstrates the influence of a fully ionized scatter- 
ing medium on the pulse shape of Her X-1 for different values 
of Tes and a variable thickness of the slab. It is obvious that for 
large optical depth values (Tes > 3), even a thin scattering layer 
((/ < 0.1) is sufficient to completely hide the pulse. In addition 
the structure of the pulse profile is steadily "washed-out" to an 
almost sinusoidal pulse shape. For Tes < 1-5 the profile's sub- 
structure (soft trailing shoulder) is still apparent even for high 
values of d. Note especially the shift in pulse phase, which is 
pronounced in the modified pulse profile for Tes > 2.0. This 
phase shift corresponds to the shift of the maximum of G(f, to) 
apparent in Fig.|9l 

5.2. Analysis of the Pulse Profile 

To simulate the variation of the pulse profile during the turn-on 
we calculated G{t, to) for the same energy ranges we used in 
Sectionl^and each RXTE orbit for two cases: 



1 . G„it, to) for a neutral corona, with NH,n 
analysis for each single orbit. 

2. Gion(f, to) for a completely ionized 
1.0 10-2 cm-2 <A?Hes < 9.5 10-4 



A^H of the spectral 



corona with 



cm 



-2 



The Green's functions Gn(f, fo) and Gion(f, to) can then be 
compared directly with the spectral model given by Eq. Q. 
Following the notation introduced in Sect. 13. II Gion(r, fo) corre- 
sponds to the model component IVIC I and Gn(f, fo) to the model 
component MC II used for the spectral analysis. For a medium 
where a fraction / is fully ionized, the total Green's function is 
given by G(f, fo) = (1 - /) Gio„(f, fo) + /Gn(f, fo). Using a fixed 
Gn(f , fo) with the respective A^e.n for each orbit from the spectral 
analysis and variable Gion(f, fo), we can simulate pulse profiles 
an observer located at infinity would see, by applying Eq. (0} 
to the "template" pulse profile. As mentioned above, the time 
scale of the simulated light curves is normahzed to the thick- 
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Fig. 11. Fit of simulated pulse profiles to observed pulse profiles for the selected observations (compare Fig.|6jl. The solid line 
represents the count rate of the observed pulse profile of the indicated orbit, and the dashed line represents the simulated best-fit 
pulse profile. The residuals of the fits are given at the bottom. All energy bands are fitted simultaneously. 



M. Kuster et al.; Her X-1 Turn-On 



11 




1.0 

Pulse Phase 

Fig. 10. Her X-1 pulse profile of orbit 30 (black) transmitted by 
a fully ionized scattering medium. The electron optical depth 
varies from top to bottom from Tes = 0.5 to Tes = 5.0. In addi- 
tion, for each optical depth value the thickness d of the scatter- 
ing medium is set to 0.02 (red), 0.02 (yellow), 0.2 (blue), and 
0.4 times (dark blue) the light crossing time of the slab. 



ness of the corona, i.e., f - fo is measured in units of the light 
crossing time of the slab. Therefore, Green's functions for dif- 
ferent values of d can be obtained by simply rescaling the time. 
For our analysis we chose d between 0. 1 and 6 light seconds, 
appropriate for the dimension of the ac cretion disk in Her X-1 
with Tin ~ 10 cm and rout ~ 10 cm dCheng et al.lll995l) . To 
obtain a proper flux normalization, the integrated flux of the 
profile of orbit 30 was set to unity. All other simulated pulse 
profiles are normalized relative to this flux. From the spectral 
fitting of the partial covering model to the observed spectra, the 
parameters of the spectral model components MC 1 and MC II 
are known. By integrating the differential photon flux of the 
model spectra over the specific energy bands of Fig.|6l the total 
flux per energy range can be calculated. Using the integrated 
photon flux both Green's functions, Gn(f, fo) and Gionit, h)), can 
then be normalized to the observed flux according to 
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Fig. 12. A^H.es as a function of time determined with the method 
of pulse profile fitting as described in the text. The error bars 
are at the Icr confidence level. 



where Nph{E) is the differential photon flux and R(h, E) the de- 
tector response matrix. Finally we performed a minimiza- 
tion fit of simulated to observed pulse profiles as shown in 
Fig-Elfor each single RXTE orbit. The different energy ranges 
were fitted simultaneously. This procedure allows us to deter- 
mine A^H.es, and d as well as their uncertainty. 

As an example, the results of the fit of simulated to ob- 
served pulses are shown in Fig.[nfor different orbits. The re- 
sult for the selected orbit demonstrates, that qualitatively and 
quantitatively the simulation can well reproduce the observed 
pulses in all energy bands. Similar results can be achieved for 
observations not shown in Fig. ^2 The integrated absolute flux 
in all energy bands is reproduced with an accuracy < 10%. For 
observations earlier than orbit 04 the analysis suffers from the 
low signal to noise ratio in the observed spectra and a resulting 
high uncertainty in the normalization of Gion(f, to) and G nit, to). 

Fig-Elshows the overaU development of the electron op- 
tical depth over the full turn on. In the figure we give Tes in 
terms of the electron column depth, A^n.es, to enable a direct 
comparison with the results of the spectral fitting. During the 
early phases of the turn on, A^e.es is very high and the absorb- 
ing medium is Compton thick. For the later observations, A^H.es 
decreases and then levels out at a constant level. These results 
are a qualitative confirmation of the behavior already expected 
from the visual inspection of the pulse profile evolution: dur- 
ing the early stages the pulse profile is completely dominated 
by scattering, while in later stages the direct radiation becomes 
more and more important. We note, however, that the overall 
electron column deduced from our Monte Carlo simulations 
is significantly higher than the column deduced from spectral 
fitting. Given the explorative character of the Monte Carlo sim- 
ulations and the simplifications introduced in the model, this 
is not unexpected. For example, our assumption of the scat- 
tering medium being a mixture of either completely neutral or 
completely ionized is certainly an oversimplification, as is the 
assumption of a slab geometry. On the other hand, as we will 
show below, despite these simplifications the overall trend of 
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A^H,es is in agreement with the common models for the 35 day 
turn-on and thus reduces the ambiguities from the spectral de- 
composition. 

6. Conclusion and discussion 

In this paper we have shown that the evolution of the X- 
ray spectrum and pulse profile during the 35 day turn-on of 
Her X- 1 can be explained by invoking a varying contribution 
of scattered and heavily absorbed photons to the observed data 
(Sect. |3and|5}- Using Monte Carlo simulations, in Sect. |5]we 
showed that the observed behavior of these components ap- 
pears to be consistent with the results of pulse profile analysis 
with theoretical Green's functions for the scattering and pho- 
toelectric absorption in the accretion disk. Despite the existing 
limitations, such as modeling the accretion disk rim by a sim- 
ple slab with uniform density or modeling the scatterer as either 
fully ionized or neutral, the methods applied here show that the 
distinct contributions to the final pulse profile and spectrum can 
in principle be separated by making use of the smearing of the 
pulse profile caused by the scattering in the plasma. Our results 
confirm earlier work based on X-ray spectral analysis alone on 
the nature of the turn-on of the 35 day cycle and on the na- 
ture of the ac cretion di sk of Her X-1 ( fPavison & Fabian' 1977^, 
|Parmaretal]il98Q: .Becker et al.lfl977: .Burwitz etal...200L) . 
Our analysis yields an optical depth of Tes ~ 3-10 for the 
scattering medium which is necessary to explain the observed 
smearing of the pulse profile. Such optical depths are consistent 
with coronal models of neutron stars in LMXRBs (e.g., iMiUeJ 
I2OOOI and references therein). 

The behavior summarized above can be explained by a sim- 
ple geometric model that takes into account the outer rim of an 
accretion disk that opens the line of sight to the neutron star and 
the influence of a hot accretion disk corona sandwiching the ac- 
cretion disk. Fig.ll3lshows the positions of the accretion disk, 
the disk corona, and the location of the observer are shown for 
difl'erent times of the turn-on. The observed PCA spectra corre- 
sponding to the phases of the turn-on and the unfolded spectral 
model components for the energy range 3 < £ < 1 2 keV are 
shown in the right panel of Fig.^] The individual components 
of the spectral model are shown separately. 

At the start of the turn-on, Thomson scattered radiation 
from the corona, which is partially absorbed, dominates the 
observed spectrum. The topmost image of Fig. fOl "Beginning 
of turn-on" illustrates the orientation of the disk and the loca- 
tion of the observer relative to the disk for this time. Photons 
following beam A, marked by a dashed line, are scattered in 
the corona and partially absorbed in the disk rim (beam C). 
The direct view to the neutron star at this time is still blocked 
and photons following beam B cannot reach the observer. The 
corresponding spectra are those observed in orbits 00 and 04. 
Since the scattered spectral model component MC I dominates 
the observed spectral flux, only a small fraction of absorbed 
flux can be detected. As a result, the total observed spectrum, 
which is the sum of all model components, shows only a weak 
signature of absorption and low values ofNn are measured. On 
the other hand, however, those scattered photons which finally 
reach the observer have undergone many scattering events such 



that the pulse profile is heavily distorted until orbit 06. In the 
pulse profile fitting, this phase of the turn-on will therefore be 
characterized by large values of A^H.es- 

As soon as the disk starts to open the line of sight to the neu- 
tron star (indicated by the downward moving disk in the image 
"Mid of Turn-On" of Fig. [131. the parts visible of the corona in- 
crease and consequently the observed flux in MC I increases as 
well. In parallel, the apparent absorbed flux increases, because 
the disk becomes more and more transparent for photons scat- 
tered in the corona (beam C). Around the mid point of the turn- 
on, the outer disk rim starts to become optically thin for pho- 
tons directly coming from the neutron star following beam D 
("Mid of turn-on"). As a consequence the flux in component 
MC II (the ratio C in Fig. |3} rises more rapidly compared to 
the flux of MC I. Since MC II determines the curvature of the 
total observable spectrum, A^h apparently increases. Since the 
pulse profile due to MC II expected to be far less smeared, pulse 
profile fitting during this episode will shown a decrease of the 
overall A^n.es- During the further evolution the contribution of 
MC II continues to increase until it dominates the spectrum. 
This event takes place right after the turning point in the track 
of A^H shown in Fig.|^. 

Shortly before the anomalous dip, the apparent A^h rapidly 
starts to decrease. Both components MC I and MC II are now 
indistinguishable (orbit 17 and later) because photoelectric ab- 
sorption and electro n scattering become negligible. Based on 
the interpretation of ICheng et al.l ( 1 19951) . this indicates that a 
larger amount of photons from the inner parts of the accretion 
disk close to the neutron star reach the observer Finally, at the 
end of the turn-on, the neutron star is directly visible when the 
main-on starts. 

In summary, models as the one outlined above seem to be 
successful in describing the overall features of the 35 day turn- 
on of Her X-1. We note, however, that our results in Sect.|5]also 
show that for very high columns, the separation into two com- 
ponents fails due to low counting statistics during early phases 
of the turn-on although the overall picture stays consistent. The 
success of the analysis presented here, however, is a first step 
towards self-consistently modeling both, the spectral evolution 
and the timing behavior of Her X- 1 . Further refinements of the 
method are still possible, e.g., when taking the absolute phase 
shift of the observed pulse profile into account and by using 
more realistic geometries. 
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Appendix A: Results of the spectral analysis 



Table A.l. Results of the spectral fitting to the Her X-1 turn-on data with all parameters set free, except the parameters where no errors are given for and the power law index 
a which was kept fix at o' = 1.068. 
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A'h: Hydrogen column density, C: Relative ratio of absorbed and scattered radiation to the unalFected radiation, Apl: Power law normalization, (photons cm^^ s^i keV^i at 1 keV), Aum'- Line 
normalization (photons cm"* s^i in the line), E^ni'- Cut-ofF energy (keV), £foid: Folding energy (keV), E^^: Gaussian line energy (keV), crp(.: Width of the Gaussian line (keV), Aune'- Line 
normalization, E^^^: Cyclotron line energy (keV), cr^^^: Width of the cyclotron line (keV), t^^^: Depth of the cyclotron line. Uncertainties are at 90 % confidence level for one interesting parameter 
(A^2 = 2.71). 
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Table A.2. Results of the spectral fitting to the Her X-1 turn-on data for the following parameters: A^h, C, Apl, and Au„e. 
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A^h: Hydrogen column density, C: Relative ratio of absorbed and scattered radiation to the unaffected radiation, Apl: Power law normalization 
(photons cm-' s-' keV-' at 1 keV), ALj„c: Line normalization (photons cm-* s-' in the line), the following parameters were fixed: the Gaussian 
emission line was fixed at 6.45 keV with a width cr of 0.45 keV, the cut off energy E^^, at 21.5 keV, the folding energy Ef^ii at 14. 1 keV, the 
power law index a at 1.068, and the cyclotron energy E^y^ at 39.4 keV with a width of 5.1 keV. Uncertainties are at 90% confidence level for 
one interesting parameter (A^^ = 2.71). 



